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Negative-index media for matter-wave optics 
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We consider the extension of optical meta-materials to matter waves and then the down scahng 
of meta-optics to nanometric wavelengths. We show that the generic property of pulsed comoving 
magnetic fields allows us to fashion the wave-number dependence of the atomic phase shift. It can 
be used to produce a transient negative group velocity of an atomic wave packet, which results into 
a negative refraction of the matter wave. Application to slow metastable argon atoms Ar*(^P2) 
shows that the device is able to operate either as an efficient beam splitter or an atomic meta-lens. 
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PACS numbers: 03.75.-b, 03.75.Be, 37.10.Gh, 42.25.-p 

Since the pioneering work of H. Lamb [l| and 
V.G. Veselago's seminal paper [2] about so-called "left- 
handed" or "meta" media for light optics, a number of 
studies have been devoted to these new media and their 
applications (negative refraction, reversed Doppler ef- 
fect, perfect lens, van der Waals atom-surface interaction, 
etc.) 0, 0, in various spectral domains [1, HIIIj some 
of them being even extended to acoustic waves Q. Such 
media are essentially characterised by a negative value of 
the optical index, which results into opposite directions 
of the wave vector k and the Poynting vector R. Our goal 
here is to extend this concept to matter waves, and the 
first arising question is the following: what should be the 
"de-Broglie optics" equivalent of those meta-materials? 
To the energy flux in electromagnetism (R vector) cor- 
responds the atomic probability flux, namely the current 
density of probability J, or equivalently the group veloc- 
ity Vg = J, where ip is the wave- function. There- 
fore, here also, one has to reverse Vg with respect to 
the wave vector k or the phase velocity. However, as 
discussed below, contrarily to what occurs in light op- 
tics where R remains directed outwards whereas k is di- 
rected towards the light source [13] , for matter waves the 
direction of the phase velocity (k) remains unchanged, 
whereas Vg is now directed towards the source. Obvi- 
ously, because of the conservation of probability, such an 
effect is necessarily a transient effect. 
In light optics, meta-media generally consist of periodic 
ensembles of micro- or nano- structures embedded in an 
ordinary material or organized into an optical band-gap 
crystal llj. Therefore their counterpart in atom optics 
is far from being obvious because atoms at mean and a 
fortiori low velocity (a few hundreds of m/s down to a 
few m/s or less) cannot penetrate dense matter. How- 
ever a possible way to act on atomic waves in vacuo is to 
use an interaction potential due to some external field, 
for instance magnetic, electric or electromagnetic fields, 
or the van der Waals field appearing outside a solid at 
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the vicinity of its surface. Indeed when a semi-classical 
description of the external atomic motion is justified, an 
inhomogeneous static potential V(r) is equivalent to an 
optical index n(r). This comes from the fact that the 
optical path accumulated along a ray C (i.e. a classi- 
cal trajectory) is given by the integral K(r)ds, where 
K{r) = k[l — V{r)/Eo]^^'^ is the local wave number, k 
and Eq being respectively the wave number and the ki- 
netic energy of the atom in absence of potential, s is the 
curvilinear abscissa along the ray p^ . This naturally 
leads us to set n{r) = [1 — y(r)/£'o]^^^- Nevertheless, in 
agreement with our previous remark about the transient 
character of the effect, such a type of potential cannot be 
a solution for our purpose because the index n is either 
real and positive, or purely imaginary in classically for- 
bidden regions, but it is nowhere real negative. This looks 
like an impasse. However our choice of potential was 
too restrictive: a much larger class of interactions is of- 
fered by position- and time-dependent potentials. Time- 
dependent potentials have already been widely used for 
devising atom optics elements in the time domain, oper- 
atin g eit her on free atom waves (modulated atomic mir- 
rors llSjl). or trapped atoms (pulsed interferometers p^ . 
dynamics of Bose-Einstein condensates [TB'l). In this let- 
ter, we show that a novel class of recently introduced po- 
tentials - i.e. comoving potentials 16] - provides us with a 
remarkably simple solution to devise negative-index me- 
dia for atomic waves. 

Comoving potentials have been previously used in several 
experiments, most of them dealing with Stern-Gerlach 
atom interferometry [T6| . They have been described in 
detail in 17[. Only the way to produce them and their 
main characteristics will be recalled here. Two identical 
planar systems of currents, periodic in space (period A) 
and symmetric with respect to x-axis (fig.l), produce in 
the vicinity of this axis a transverse field, e.g. parallel 
to y, periodic in x with the same period. Actually only 
the lowest spatial frequency will have a significant effect. 
Then the field can be assumed to be proportional to, 
for instance, cos(27ra;/A). Now if the circuit is supplied 
with an A.C. intensity of frequency v, then the result- 
ing field is proportional to cos(27ri/(t — tg)) cos(27ra;/A) = 
i [cos{2n{L't -to)-x /A) + cos{2Tr{v{t-to) + x /A)] , where 
to is some reference time. Only the first term propagates 
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FIG. 1: Scheme of the device generating a transverse mag- 
netic field "comoving" along the x axis. It consists of two 
identical planar periodic circuits (spatial period A). When 
the circuits are supplied with an A.C. current of frequency v 
, they produce along x a field that propagates at the velocity 
u = vK. When a frequency spectrum H(t/ ) is used, a pulsed 
field is generated (see text). 



in the same direction as the atoms {z > 0), at velocity 
u = vK. The period A being fixed, the velocity u can 
be varied via v. In particular it can be tuned on the 
atomic velocity, giving rise to an especially important 
phase-shift on the atomic wave contrarily to the second 
counter-propagating term which can be ignored. When a 
continuous spectrum ii(y) is used instead of a single fre- 
quency, the resulting (pulsed) interaction potential takes 
the general factorized form: 



V[t,x) = s{t) cos (Stt^) 



(1) 



The time-dependent factor is real and can be expressed 
as: s(t) = g \XBMBf{t). Here g is the atomic Lande fac- 
tor, [juB the Bohr magneton, M the magnetic quantum 
number and B the maximum of the field magnitude, /(t) 
being a reduced dimensionless signal the maximum am- 
plitude of which is 1. Since s(t) is real, its frequency 
spectrum Hiy) obeys the relation: 



(2) 



where * means complex conjugate. The motion along the 
X axis of a wave packet submitted to the potential V[t, x) 
is governed by the time-dependent Schrodinger equation: 



ih dtip — 



2m 



(3) 



where m is the atomic mass. Taking the Fourier trans- 
form of eq. ([3]) for variables x, k {k being the wave num- 
ber) one gets: 



ihdtC{t,k) 
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where C and W are the spatial Fourier transforms of V' 
and V, and ® the convolution product. Setting C{t, k) = 
exp(— i^t)r(t, fc), one readily gets: 
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FIG. 2: Trajectory of the wave packet centre (coordinates 
^(t), C,(t) along X- and 2:-axes) in the comoving magnetic field 
(dashed line) propagating in the iix direction in the region 
z > Q. The entrance plane is perpendicular to the z-axis. 
Rays start from a point-like source 5* located at 2 cm from 
the entrance plane. These rays are in plane x, z. They ini- 
tially make with the z axis different incident angles ranging 
from to 0.12 rad. All rays exhibit a negative refraction and 
finally emerge, for z > 24 mm, parallel to their initial direc- 
tion. Calculations are made for Ar*(^P2, M=2) metastable 
atoms, the velocity of which is wo = 20 m/s. The maximum 
magnitude of the magnetic field (see text) is 400 Gauss. Dis- 
tances are in meters. Note the difference in scale for x and z 
axis. 



where k — 2tt/ A {k << k). A simplification arises be- 
cause the k-derivatives of F are largely dominated by 
those of the exponential factor exp(— i|^t) as soon as 
the wave packet has moved over a distance large com- 
pared to its own width [T7| . Under such conditions: 

ih dtVit, k) « s{t) cos ( 27r^i j r(t, k) (6) 



Then: 



r(t, fc) « r(0, fc) exp[i(/3(fc,t)] 



(7) 



where the phase shift Lp is given by : 



^(t, k) = -ft-i ^ dt' s{t') cos (2Tr^t'^ (8) 

This real phase-shift takes a limiting value at t infinite, 
namely, assuming a perfect synchronisation of the wave 
packet with the field pulse (eqs. 1 and 2): 



hk ^ 



(9) 



This can be seen as a genericity property of comoving 
fields, in the sense that, in principle, any k-dependence of 
the phase-shift can be fashioned, using a convenient H{v) 
spectrum. As suggested previously in it could be 
used to balance the natural spreading of a wave packet. 
The factor s{t), or its frequency spectrum H{i>), being 
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given, one may derive the evolution of the wave packet 
in the corresponding potential. In particular, the semi- 
classical motion of the wave packet center is derived from 
the stationary phase condition: 

hk^ hk 

dk [kx - — t + ip{k,t)]^x t + dkfik,t) ^0 (10) 

2m m 

Then the abscissa of the center as a function of time is 
simply: 



m 



(11) 



hkf) being the central momentum value. The quantity 
(5^(t) — —dkf(k,t)\kg represents the spatial shift induced 
by the potential. It is given by: 



S^(t) = /* dt' t'sit') sin ( 2Tr-^A (12 

mA Jq \ niA J 



At large values of t, when the potential pulse is over, this 
shift naturally tends to a definite limiting value, namely: 



mA mA 



(13) 



where H' is the derivative of H. The group velocity along 
X is readily derived from (11) and (12): 



Vgx {t,ko) 



— - dtdk(p{k,t)\ko 
m 

hkn 27r , , / hkn , 
T:^-3T*^Wsin 2^-|i (14) 
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The quantity f q — hk^ j (mA) can be called the "reso- 
nance" frequency in the sense that the velocity u{i'o) = 
vqA of this spectral component of the comoving field co- 
incides with the atomic group velocity. For i < as 
well as for large t values, the wave packet recovers the 
group velocity hko/m of the free propagation. In be- 
tween, it is clearly seen from (14) that, by a proper choice 
of s(i), negative group velocities can be (transiently) ob- 
tained whereas the phase velocity remains positive. Such 
behaviour is the signature of a left-handed or "meta"- 
medium for atomic waves. It gives rise to new interest- 
ing effects. Let us assume that the entrance interface for 
atoms is the Oxy plane. Then the trajectory related to 
a positive magnetic quantum number (M > 0) is sub- 
jected to a negative refraction, another signature of a 
left-handed or "meta" medium, whereas for M < the 
trajectory undergoes an ordinary refraction, with an ef- 
fective index smaller than 1. In the present case, the 
medium is anisotropic. It is uni-axial and optically ac- 
tive. Fig. 2 shows rays issued from a point-like source S 
located at 2 cm from the interface, experiencing a neg- 
ative refraction (M = -1-2) in the xz plane, for different 
incidence angles 9 ranging from to 0.12 rad. Numerical 
calculations have been carried out in the case of a nozzle 
beam of metastable argon atoms Ar*(^P2), slowed down 
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FIG. 3: Position of the image S' , where the support of the 
emerging ray crosses the z axis. The distance AZ = SS' is 
shown (in meters) as a function of the incident angle (in rad) . 
At no = 20 m/s, the stigmatism is realized with an accuracy 
better than 90% for incident angles lower than 0.04 rad. The 
chromatic dispersion is also shown (lower curve : no = 18 
m/s, mid curve : wo = 20 m/s, upper curve : «o = 22 m/s). 
Note the absence of chromatism at an incidence angle of 0.09 
rad. 




FIG. 4: (a) 3D representation of half a cone of rays issued from 
a point-like source and making an angle 9 = 0.1 rad with the z 
axis. All rays exhibit a negative refraction and finally emerge 
parallel to their initial direction, (b) Same as (a), with a 
2D comoving potential s{t) = [cos(27ra;/A) -I- cos(27ry/A)] (see 
text). The system behaves as a spherical meta-lens. Note that 
this surface slightly differs from that generated by the curve 
^{z) rotated around the z-axis. 



at a velocity vq = 20 m/s by means of a Zeeman slower 
[3[l^. The parameters used in these calculations are, 
for the magnetic potential, A = 5 mm, B = 400 Gauss, 
f{t) = e^{t + e)-2e-*/^ for < t < n, = elsewhere, 
with e = 7.4 ms, t — 0.37 ms, ti = 1.2 ms. Actually the 
negative refraction only appears when the magnitude of 
B exceeds a threshold value (280 Gauss in the present 
case) . Rays symmetric of the previous ones with respect 
to the z axis experience opposite refractions due to the 
opposite component of the comoving field. The device 
thus behaves as a parallel plate, or a cylindrical meta- 
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FIG. 5: Effect of three subsequent pulses s{t) + 1.5s(t~lms) + 
s{t — 2ms) of comoving potential (shown in the inset), on the 
trajectories of the wave packet centre at incidence angles 6 
ranging from to 16 mrad. A refocusing is seen, with a final 
transverse astigmatism lower than 25/im. 

lens which gives an image S' of the source point S. In a 
real experimental situation, e.g. in a cooled thermal or 
nozzle beam, there exists an atomic velocity distribution. 
The angular velocity spreading is related to the stigma- 
tism of the system, as shown in fig. 3. For an angular 
aperture of 0.08 rad, the axial stigmatism is limited to 
3.3 mm. The dispersion of the velocity modulus gives rise 
to a chromatic effect. This effect is also shown in fig. 3, for 
a wide variation (±10%) of uq. It is maximum (±4 mm) 
at 6* = and cancels sX 9 ~ 0.09 rad. Such "thermal" 
spreadings are much larger than those experimentally ac- 
cessible {59 < 0.01 rad, \5vo/vo\ « 1 - 2%). Therefore 
effects related to negative group velocities should be ob- 
servable. Obviously ultra-cold atoms, and a fortiori, con- 
densates are expected to give rise to even more striking 
effects. The evolution of a conical map of rays with — 
0.1 rad is shown in fig. 4a. The previous treatment can 
be easily generalized to a 2D potential, provided it allows 
a separation of variables x and y, such as a potential in 
cos(27ra:/A) -|-cos(27r2//A). Under such conditions the de- 



vice behaves at small incidence angles as a spherical lens 
(see fig. 4b). Atom wave re- focusing is actually predicted 
when using a series of two or three comoving pulses (Fig. 
5). It is also worth noting that, with atoms initially spin- 
polarized in a linear superposition of Zceman sub- levels, 
ID-comoving potentials could be used in atom interfer- 
ometry, as very efficient beam splitters. 

Finally, along the similarities between meta-media in 
light optics and atom optics, we should underline some 
of their basic differences. They originate in the distinct 
characteristics of wave propagation for Maxwell waves 
and "de Broglie - Schrodinger" waves: contrary to elec- 
tromagnetic waves, de Broglie waves for massive particles 
undergo a longitudinal wave packet spreading due to the 
fundamental vacuum dispersion in matter optics (as im- 
posed by the dispersion relation lo = hk"^ /2m). Contrar- 
ily to optical meta-materials, in scalar de Broglie optics, 
atomic meta-media with negative index are characterised 
by a reversal of the group velocity, the direction of the 
wave vector remaining unchanged. In general, retarda- 
tion effects in meta-optics are important. On the other 
hand, in matter optics, the response time is mainly gov- 
erned by the finite velocity of the massive atomic wave. 
In conclusion, atom optics with so-called meta-media 
paves the way to a new field of original operations to 
be performed on matter waves [20] . If the negative index 
for matter waves can be tuned close to —1, then 100% 
transmission is expected at normal incidence (zero Fres- 
nel reflection coefficient). In analogy to light wave optics 
I3:] , one might be able to devise atomic lenses allowing one 
to re-focus atom waves without diffraction limitations. 
The latter characteristics should appear provided that 
evanescent matter waves can be properly reconstructed 
inside the atomic meta-medium. This key property for 
sub "de-Broglie-wavelength" focussing will be considered 
in a forthcoming publication. This would open novel ap- 
plications in atom nano-lithography. 
Authors are members of the Institut Francilien de 
Recherche sur les Atomes Froids (IFRAF). 
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